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This Letter reports a measurement of the flux and energy spectrum of electron antineutrinos from six 
2.9 GW th nuclear reactors with six detectors deployed in two near (effective baselines 512 m and 561 m) and 
one far (1,579 m) underground experimental halls in the Daya Bay experiment. Using 217 days of data, 296,721 
and 41,589 inverse beta decay (IBD) candidates were detected in the near and far halls, respectively. The mea¬ 
sured IBD yield is (1.55 ± 0.04) x 10” 18 cm 2 /GW/day or (5.92 ± 0.14) x 10 -43 cm 2 /fission. This flux 
measurement is consistent with previous short-baseline reactor antineutrino experiments and is 0.946 ± 0.022 
(0.991 ± 0.023) relative to the flux predicted with the Huber+Mueller (ILL+Vogel) fissile antineutrino model. 
The measured IBD positron energy spectrum deviates from both spectral predictions by more than 2 o over the 
full energy range with a local significance of up to ~4 o between 4-6 MeV. A reactor antineutrino spectrum of 
IBD reactions is extracted from the measured positron energy spectrum for model-independent predictions. 

PACS numbers: 14.60.Pq, 29.40.Mc, 28.50.Hw, 13.15.H-g 
Keywords: antineutrino flux, energy spectrum, reactor, Daya Bay 


Reactor antineutrino experiments have played a key role in 
developing the picture of neutrinos in the Standard Model of 
particle physics. They provided the first experimental obser¬ 
vation of (anti)neutrinos Q], confirmed neutrino oscillation as 
the solution to the solar neutrino problem 12, provided the 
first detection of geo-neutrinos 0, and discovered the non¬ 
zero neutrino mixing angle () \:>, |4][5]. Forthcoming reactor an¬ 
tineutrino experiments are aiming to further explore the nature 
of neutrinos by determining the neutrino mass hierarchy, pre¬ 
cisely measuring neutrino mixing parameters, and searching 
for short-baseline neutrino oscillation 0. Over the last five 
decades, reactor antineutrino experiments have measured the 
flux and spectrum of antineutrinos at various distances from 
nuclear reactors ranging from ~10 m to several hundred kilo¬ 
meters. These measurements were found to be in good agree¬ 
ment 0 with predictions derived from the measurements of 
the beta spectra at ILL 181 4101 and Vogel’s theoretical calcu¬ 
lation mi when considering the effect of three-neutrino os¬ 
cillation. In 2011, re-evaluations of the reactor antineutrino 
flux and spectrum with improved theoretical treatments were 
carried out OH S3, and determined the flux to be higher than 
the experimental data. This discrepancy is commonly referred 
to as the “Reactor Antineutrino Anomaly’ ’Da and may be a 
sign of new physics or insufficient fissile antineutrino model¬ 
ing. Precision measurements by modern reactor antineutrino 
experiments can shed light on this issue and probe the physics 
underlying current reactor antineutrino predictions. An ac¬ 
curate determination of the reactor antineutrino spectrum can 
also provide valuable input to next-generation single-detector 
reactor antineutrino experiments D3- 


This Letter reports measurements of the reactor antineu¬ 
trino flux and spectrum based on 217 days of data from the 
Daya Bay experiment. The Daya Bay reactor complex con¬ 
sists of three nuclear power plants (NPPs), each hosting two 
pressurized-water reactors. The maximum thermal power of 
each reactor is 2.9 GW th- The data used for this analysis 
comprises 338,310 antineutrino candidate events collected in 
six antineutrino detectors (ADs) in the two near experimental 
halls (effective baselines 512 m and 561 m) and the one far 
hall (effective baseline 1,579 m). This is the largest sample of 
reactor antineutrinos, comparable to that from the BUGEY-4 
experiment Da. A more detailed description of the experi¬ 
mental setup and the data set is given in Ref. Da- 

In reactor cores, electron antineutrinos (v e ) are emitted 
isotropically from fission products of four primary isotopes: 
235 U, 238 U, 239 Pu, and 241 Pu. The number of v e with energy 
E emitted from a reactor at a time t can be predicted using 


d 2 (f>(E, t) 
dEdt 


E fi(t)Si(E)c^{E,t) + SsME,t), 

G) 


where the sums are over the four primary isotopes, W t h(t) 
is the reactor thermal power, /j(t) is the fraction of fis¬ 
sions due to isotope i, e, is the average thermal energy re¬ 
leased per fission, Si(E) is the v e energy spectrum per fis¬ 
sion, c" e (i?,t) is the correction to the energy spectrum due 
to reactor non-equilibrium effects of long-lived fission frag¬ 
ments, and SsNF(E,t) is the contribution from spent nuclear 
fuel (SNF). At Daya Bay, the NPPs monitor the reactor power 
in real-time and simulate the evolution of the fuel composi- 
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tion using the SCIENCE software package mm. The mea¬ 
sured power (0.5% uncertainty USES) and simulated fis¬ 
sion fractions (~5% relative uncertainty ll23l ) of each core 
are provided to the Daya Bay collaboration. Simulation of re¬ 
actor cores based on DRAGON ll24l was constructed to study 
the correlations among the fission fractions of the four iso¬ 
topes ll25ll . The energies released per fission (0.2-0.5% un¬ 
certainty) were from Ref. (26) . Non-equilibrium (30% uncer¬ 
tainty) and SNF (100% uncertainty) corrections were applied 
following Refs. JT2] and (23 (28), with -0.5% and -0.3% 
contributions to the total antineutrino rate, respectively. Com¬ 
bining the uncertainties of reactor power, fission fractions, 
and non-equilibrium and SNF corrections, the total reactor- 
uncorrelated uncertainty of antineutrino flux is 0.9%. 

Two fissile antineutrino spectrum models were used for 
Si(E) in Eq. (JTJ) to predict the reactor antineutrino flux and 
spectrum. The ILL+Vogel model refers to the conventional 
ILL model ifSTTfOl of 235 U, 239 Pu, and 241 Pu, and the theo¬ 
retical model of 238 U from Vogel fill . The Huber+Mueller 
model refers to the recent re-evaluation of 235 U, 239 Pu and 
241 Pu from Huber lfl3l . and that of 238 U from Mueller et 
al. fl2l . The Huber+Mueller model was chosen as a refer¬ 
ence because of its improved theoretical treatments in beta- 
to-antineutrino conversions, and the information it provides 
about uncertainties and their correlations. 

Reactor antineutrinos were detected via inverse beta decay 
(IBD) reactions in the gadolinium-doped liquid scintillator 
(GdLS) of the Daya Bay ADs. The total number of detected 
IBD events T in a given AD was estimated as 


t = E 


Np £ibd 


JjP SUI (E, Li)a 1 B n(E) ^ dEdt, 

( 2 ) 


where d 2 <pi(E,t)/dEdt is the differential antineutrino rate 
from the i- th reactor core given in Eq. Q, CTIBD {E) is the 
cross section of the IBD reaction, Lj is the distance between 
the center of the detector and the i-th core {29), P S \ir{E, Li) is 
the survival probability due to neutrino oscillation, Np is the 
number of target protons {30), and £ibd is the efficiency of 
detecting IBD reactions. The cross section <tibt>{E) was eval¬ 
uated based on the formalism in Ref. ED. Physical constants 
including the neutron lifetime (880.3 ± 1.1 s) were taken from 
the Particle Data Group (32). 

IBD candidates were selected by requiring a time coinci¬ 
dence between a prompt signal from an IBD positron includ¬ 
ing its annihilation energy, and a delayed signal from an IBD 
neutron after capturing on Gd, as described in Refs. mm 
The energy of interacting antineutrinos, E, is closely related 
to the prompt energy of the IBD positrons, /i prolnpt : without 
detector effects, E pmmpt ~ E + ( M p — M n — M e ) + 2M e = 
E — 0.78 MeV, where M p , M n , and M e are the proton, 
neutron, and electron masses. The reported dataset includes 
296,721 and 41,589 IBD candidates at the near and far halls, 
respectively. Corresponding background rates and spectra 
were estimated in Ref. fl7) . with about 5,470 ± 240 and 1,894 
± 43 background events at the near and far halls, respectively. 


TABLE I. Summary of IBD selection efficiencies and their AD- 
correlated uncertainties. The uncertainties are given in relative units. 



Efficiency (e) 

Uncertainty (<5 Je) 

Target protons 

- 

0.47% 

Flasher cut 

99.98% 

0.01% 

Capture-time cut 

98.70% 

0.12% 

Prompt-energy cut 

99.81% 

0.10% 

Gd-capture fraction 

84.17% 

0.95% 

Delayed-energy cut 

92.71% 

0.97% 

Spill-in correction 

104.86% 

1.50% 

Combined 

80.6% 

2.1% 


The relative analysis of IBD rates for sin 2 29±3 requires 
estimates of uncertainties that are uncorrelated among ADs 
while the measurement of flux is dominated by uncertainties 
that are correlated among ADs. A detailed study of the event 
selection efficiencies was carried out using Monte Carlo sim¬ 
ulation (MC)-data comparisons with the Daya Bay simulation 
framework based on Geant4 (33) . The previous study of 
efficiencies is described in detail in Ref. {29). Estimates of 
each selection efficiency and detector characteristic are sum¬ 
marized in Tableland are briefly described in order below. Ef¬ 
ficiencies of flashing-photomultiplier tube (flasher), capture¬ 
time, and prompt-energy selections were determined as de¬ 
scribed in Ref. l29l utilizing an updated Daya Bay IBD MC. 
The IBD neutron-Gd capture fraction is dependent on the tar¬ 
get’s Gd concentration and on the escape, or ‘spill-out’, of 
IBD neutrons from the target. The former has been measured 
using neutron calibration sources deployed at the detector cen¬ 
ter, while the latter was estimated with MC-data comparisons 
of source deployments throughout the GdLS volume with a 
manual calibration system (34). The efficiency for detecting 
Gd-capture IBD neutrons, also called the delayed energy cut 
efficiency, is dependent on the amount of Gd-capture 7 energy 
deposited outside the scintillator, and was determined using 
MC benchmarked to the IBD Gd-capture spectrum from data. 
Finally, in order to account for contributions from IBD in¬ 
teractions outside the GdLS target, we have applied a spill-in 
correction determined using MC-data comparisons of IBD co¬ 
incidence time and reconstructed position distributions. The 
updated detector efficiency e was estimated to be 80.6% with 
an AD-correlated fractional uncertainty 8Je of 2.1%. Appli¬ 
cation of additional AD-dependent muon veto and multiplicity 
cut efficiencies, described in detail in Refs. {HIED, produced 
total detection efficiencies £ibd ranging from 64.6% to 77.2% 
among ADs. The total correlated uncertainty was dominated 
by the spill-in correction, whose uncertainty enveloped the in¬ 
dividual uncertainties provided by three independent methods 
and was limited by small biases in position reconstruction. A 
cross-check of the spill-in effect provided by data-MC com¬ 
parisons of neutron sources deployed outside the target vol¬ 
ume showed agreement well within this uncertainty. 

To extract the rate of IBD interactions at Daya Bay, the 6 13 - 
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driven oscillation effect must be corrected for in each detec¬ 
tor. A normalization factor R was defined to scale the mea¬ 
sured rate to that predicted with a fissile antineutrino spectrum 
model. The value of R, together with the value of sin 2 20 13 , 
were simultaneously determined with a \ 2 similar to the one 
used in Ref. EQ: 

2 _ V - ' [Md R ’ Td( 1 ± £d ± Sr + e d) + VdY 
X M d + B d 



where M d is the number of measured IBD events in the d-th 
detector with backgrounds subtracted, B d is the correspond¬ 
ing number of background events, T d is the number of IBD 
events predicted with a fissile antineutrino spectrum model 
via Eq. (|2j, and ujf is the fractional IBD contribution from 
the r-th reactor to the d-th detector determined with baselines 
and reactor antineutrino rates, ay (0.9%) is the uncorrelated 
reactor uncertainty, a d (0.2% El) is the uncorrelated de¬ 
tection uncertainty, ag d is the background uncertainty listed 
in Ref. 07], and (td (2.1%) is the correlated detection un¬ 
certainty, i.e. the uncertainty of detection efficiency in Ta¬ 
ble [I] Their corresponding nuisance parameters are ay, e d , 
r)d, and ey>, respectively. The best-fit value of sin 2 2 di 3 = 
0.090 ± 0.009 is insensitive to the choice of model. The best- 
fit value of R is 0.946±0.022 (0.991 ±0.023) when predicting 
with the Huber+Mueller (ILL+Vogel) model. Replacing the 
Mueller 238 U spectrum with the recently-measured spectrum 
in Ref. Il35l yields negligible change in R. The uncertainty in 
R is dominated by the correlated detection uncertainty a d. 

With the oscillation effect for each AD corrected using 
the best-fit value of sin 2 2 ±13 in Eq. ([ 3 J, the measured IBD 
yield for each AD is expressed in two ways: the yield per 
GW th per day, Y, and equivalently, the yield per nuclear fis¬ 
sion, af. These results are shown in the top panel of Fig. |T] 
The measured IBD yields are consistent among all ADs after 
further correcting for the small variations of fission fractions 
among the different sites. The average IBD yield in the three 
near ADs is Y = (1.55 ± 0.04) x 10 -18 cm 2 /GW/day, or 
07 = (5.92 ± 0.14) x 10 -43 cm 2 /fission. These results are 
summarized in Table [H| along with the flux-weighted average 
fission fractions in the three near ADs. 

A global fit for R was performed to compare with previous 
reactor antineutrino flux measurements following the method 
described in Ref. |(36l . Nineteen past short-baseline (<100 m) 
measurements were included using the data from Ref. 031 . 
The measurements from CHOOZ l(37l and Palo Verde l38l 
were also included after correcting for the effect of standard 
three-neutrino oscillations. All measurements were compared 
to the Huber+Mueller model. All predictions were fixed at 
their nominal value in the fit. The resulting past global average 
is R p g ast = 0.942±0.009 (exp.)±0.025 (model). Daya Bay’s 
measurement of the reactor antineutrino flux is consistent with 
the past experiments. Including Daya Bay in the global fit, the 


TABLE II. Average IBD yields (Y and 07 ) of the near halls, flux nor¬ 
malization with respect to different fissile antineutrino model predic¬ 
tions, and flux-weighted average fission fractions of the near halls. 


IBD Yield 

Y ( cm 2 /GW/day) 

(1.55 ±0.04) x 10 " 18 

<jf (cm 2 /fission) 

(5.92 ±0.14) x 10 -43 

Data / Prediction 

R (Huber+Mueller) 

0.946 ± 0.022 

R (ILL+Vogel) 

0.991 ± 0.023 

235|j . 238 -q . 239. 241p u 

0.586 : 0.076 : 0.288 : 0.050 


new average is R g = 0.943 ± 0.008 (exp.) ± 0.025 (model). 
The results of the global fit are shown in the bottom panel of 

Fig- |U 

Extending the study from the integrated flux to the en¬ 
ergy spectmm, the measured prompt-energy spectra of the 
three near-site ADs were combined after background subtrac- 



Distance (m) 


FIG. 1. Top: Rate of reactor antineutrino candidate events in the six 
ADs with corrections for 3-flavor oscillations (closed circles), and 
additionally for the variation of flux-weighted fission fractions at the 
different sites (open squares). The average of the three near detectors 
is shown as a gray line (and extended through the three far detectors 
as a dotted gray line) with its la systematic uncertainty (gray band). 
The rate predicted with the Huber+Mueller (ILL+Vogel) model and 
its uncertainty are shown in blue (orange). Bottom: The measured 
reactor P e rate as a function of the distance from the reactor, nor¬ 
malized to the theoretical prediction with the Huber+Mueller model. 
The rate is corrected for 3-flavor neutrino oscillations at each base¬ 
line. The blue shaded region represents the global average and its la 
uncertainty. The 2.7% model uncertainty is shown as a band around 
unity. Measurements at the same baseline are combined for clarity. 
The Daya Bay measurement is shown at the flux-weighted baseline 
(573 m) of the two near halls. 
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tion and compared with predictions. The antineutrino spec¬ 
trum at each detector was predicted by the procedure de¬ 
scribed above, taking into account neutrino oscillation with 
sin 2 20 13 = 0.090 and Ara 2 e = 2.59 x 10 -3 eV 2 based on 
the oscillation analysis of the same data El. The detector 
response was determined in two ways. The first method se¬ 
quentially applied a simulation of energy loss in the inactive 
acrylic vessels, and analytical models of energy scale and en¬ 
ergy resolution. The energy scale model was based on empir¬ 
ical characterization of the spatial non-uniformity and the en¬ 
ergy non-linearity with improved calibration of the scintillator 
light yield and the electronics response l39l . The uncertainty 
of the energy scale was about 1% in the energy range of reac¬ 
tor antineutrinos l39l . The second method used full-detector 
simulation in which the detector response was tuned with the 
calibration data. Both methods produced consistent predic¬ 
tions for prompt energies above 1.25 MeV. Around 1 MeV, 
there was a slight discrepancy due to different treatments of 
IBD positrons that interact with the inner acrylic vessels. Ad¬ 
ditional uncertainty below 1.25 MeV was included to cover 
this discrepancy. 

Figure [2] shows the observed prompt-energy spectrum and 
its comparison with the predictions. The spectral uncertainty 
of the measurement is composed of the statistical, detector 
response and background uncertainties. Between 1.5 and 7 
MeV, it ranges from 1.0% at 3.5 MeV to 6.7% at 7 MeV, and 
above 7 MeV it is larger than 10%. The predicted spectra were 
normalized to the measurement thus removing the dependence 
on the total rate. Agreement between a prediction and the data 
was quantified with the y 2 defined as 


x 2 = ^(JV“ bs - Nf red ) V," 1 (N° hs - Af ed ), (4) 

i,3 
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FIG. 2. Top panel: Predicted and measured prompt-energy spectra. 
The prediction is based on the Huber+Mueller model and normalized 
to the number of measured events. The highest energy bin contains 
all events above 7 MeV. The gray hatched and red filled bands rep¬ 
resent the square-root of diagonal elements of the covariance matrix 
WVu) for the reactor related and the full (reactor, detector and back¬ 
ground) systematic uncertainties, respectively. The error bars on the 
data points represent the statistical uncertainty. Middle panel: Ra¬ 
tio of the measured prompt-energy spectrum to the predicted spec¬ 
trum (Huber+Mueller model). The blue curve shows the ratio of the 
prediction based on the ILL+Vogel model to that based on the Hu¬ 
ber+Mueller model. Bottom panel: The defined y 2 distribution (y,) 
of each bin (black dashed curve) and local p-values for 1-MeV en¬ 
ergy windows (magenta solid curve). See the text for the definitions 
of these quantities. 


where _/y° bs ( prcc 9 j s t b e observed (predicted) number of events 
at the *-th prompt-energy bin and V is the covariance ma¬ 
trix that includes all statistical and systematic uncertainties. 
The systematic uncertainty portion of the covariance matrix 
V was estimated using simulated data sets with randomly 
fluctuated detector response, background contributions, and 
reactor-related uncertainties, while the statistical uncertainty 
portion was calculated analytically. A comparison to the Hu¬ 
ber+Mueller model yielded a y 2 /NDF, where NDF is the num¬ 
ber of degrees of freedom, of 43.2/24 in the full energy range 
from 0.7 to 12 MeV, corresponding to a 2.6 ct discrepancy. The 
ILL+Vogel model showed a similar level of discrepancy from 
the data at 2.4 ct. 

The ratio of the measured to predicted prompt-energy spec¬ 
tra is shown in the middle panel of Fig. [2] A discrepancy 
is apparent around 5 MeV. Two approaches were adopted to 
evaluate the significance of local discrepancies. The first was 
based on the y 2 contribution of each energy bin, which is eval¬ 


uated by 


jyobs _ y P re< t 


Xi = 


jyobs _ y pred 




(5) 


where y 2 , = (N° hs - N^V^ 1 (N° bs - iV pred ). 


As shown in the bottom panel of Fig. [2] there is a larger contri¬ 
bution around 5 MeV. In the second approach, the significance 
of deviations are conveyed with p-values calculated within lo¬ 
cal energy windows. A free-floating nuisance parameter for 
the normalization of each bin within a chosen energy window 
was introduced to the fitter that was used in the neutrino os¬ 
cillation analysis. The difference in the minimum y 2 before 
and after introducing these nuisance parameters was used to 
evaluate the p-value of the deviation from the theoretical pre¬ 
diction within each window. The p-values within 1-MeV en¬ 
ergy windows are shown in the bottom panel of Fig. [2] The 
p-value for a 2-MeV window between 4 and 6 MeV reached a 
similar minimum of 5.4 x 10~ 5 , which corresponds to a 4 .Oct 
deviation. The ILL+Vogel model showed a similar level of 
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discrepancy between 4 and 6 MeV. 

The number of events in excess of the predictions in the 4- 
6 MeV region was estimated to comprise approximately 1% 
of all events in both the near and far detectors. This excess 
is approximately 10% of events within the 4-6 MeV region. 
This discrepancy was found to be time-independent and cor¬ 
related with reactor power, therefore disfavoring hypotheses 
involving detector response and unknown backgrounds. A re¬ 
cent ab-initio calculation of the antineutrino spectrum showed 
a similar deviation from previous predictions in the 4-6 MeV 
region, and identified prominent fission daughter isotopes as 
a potential explanation l40l . A number of tentative explana¬ 
tions based on the nuclear physics of beta decays and fission 
yields have been put forward and are under active investiga¬ 
tion; for examples, see Refs. I40H441 . These studies suggest 
an increased uncertainty in both the yields and spectra of the 
fissile antineutrino models, which may also account for the 
discrepancy. 

From the measured IBD prompt spectrum at Daya Bay, we 
have obtained a reactor antineutrino spectrum of IBD reac¬ 
tions that can be used to make model-independent predictions 
of reactor antineutrino flux and spectra P31 . The spectrum 
was obtained by first summing the prompt-energy spectra of 
the three near site ADs weighted with their target mass rel¬ 
ative to the average target mass of all near-site ADs, M: 

^combined(-^prompt) — 5) ( ^prompt ) TF/A/, . Detector le- 

sponse effects were then removed by unfolding the combined 
prompt spectrum ^combined(fJprompt) to an antineutrino spec¬ 
trum of IBD reactions, S'combined(-E'). Finally, oscillation ef¬ 
fects were removed and each bin of the antineutrino spectrum 
was normalized to cm 2 /fission/MeV using the thermal power 
Wthit ) and fission fraction /*(f) information of each core. 
The reactor antineutrino spectrum is expressed as 


S K 


r(E) = T? 


s< 


combined 


(E) 


Psur(E) ■ Np • Ft ota i ’ 


( 6 ) 


where P ml (E) is the flux-weighted average of the survival 
probabilities P sur (E, Lt^) from the six reactors ( i ) to the 
three detectors ( d ), Np is the number of target protons in 
M, and F tota i is the total number of fissions from the sum of 
the fissions of the six reactors to the three detectors weighted 
with £iBD,d/47rL 2 d . Correcting the unfolded spectrum with 
an average survival probability resulted in a negligible bias 
(< 0.01 %). 

Detector response effects were removed with the Singular 
Value Decomposition (SVD) unfolding method ll46l . Statis¬ 
tical and systematical uncertainties are naturally propagated 
in the SVD method. The bias of unfolding was estimated by 
using detector response matrices constructed from the two dif¬ 
ferent detector response models and by using a variety of in¬ 
put antineutrino spectra which covered the uncertainties of the 
two models and those estimated in Ref. ESI. The bin-to-bin 
bias between 2.2 and 6.5 MeV was about 0.5%, which was a 
few times smaller than the statistical uncertainty. The bias out¬ 
side this region was about 4% and increased with energy due 


0.2 



Antineutrino Energy (MeV) 


FIG. 3. Top panel: The extracted reactor antineutrino spectrum and 
its correlation matrix. Bottom panel: Ratio of the extracted reac¬ 
tor antineutrino spectrum to the Huber+Mueller prediction. The er¬ 
ror bars of the data points are the square-roots of the diagonal ele¬ 
ments of the data covariance matrix, which included statistical and 
systematic uncertainties, as well as bias and the AD-correlated un¬ 
certainty from Table|l] The solid red band represents the square-roots 
of the diagonal elements of the prediction covariance matrix, which 
included reactor and Huber+Mueller model uncertainties. The hori¬ 
zontal dashed line represents the normalization factor R = 0.946. The 
vertical dashed line denotes that above 8 MeV, the Huber+Mueller 
model was extrapolated. 


to the decrease of events. The bias values were assigned as 
additional uncertainties to the unfolded spectrum. Unfolding 
performed with the Bayesian iteration method |47] 148 1 pro¬ 
duced consistent results. Between 2 and 7.5 MeV, the spec¬ 
tral uncertainty of the unfolded spectrum S'combined(-E’) ranges 
from 1.1% at 4.25 MeV to 9.3% at 7.5 MeV, and around 10 
MeV is more than 20% due to low statistics. The obtained 
reactor antineutrino spectrum and its correlation matrix are 
shown in the top panel of Fig. [3] Between 2 and 7.5 MeV, 
the uncertainty of the diagonal elements ranges from 2.7% at 
4.25 MeV to 10.4% at 7.5 MeV. The bottom panel of Fig. [3] 
is the ratio of the extracted reactor antineutrino spectrum to 
the prediction using the fissile antineutrino spectra of the Hu¬ 
ber+Mueller model and the average fission fractions listed in 
Table [B] The integral of the ratio is equal to the flux normal¬ 
ization factor R given in Table[U] The integral of the spectrum 
is equal to the yield cry given in Table [ll) The discrepancy be¬ 
tween 5 and 7 MeV corresponds to the discrepancy between 4 
and 6 MeV in the IBD prompt-energy distribution in Fig. [2] 

In summary, the Daya Bay experiment collected more 
than 330,000 antineutrino events in the data-taking period 
with six antineutrino detectors. The measured IBD yield is 
(1.55 ± 0.04) x 10~ 18 cm 2 /GW/day or (5.92 ± 0.14) x 
1 O' 13 cm 2 /fission. This flux measurement is consistent with 
the global average of previous short baseline experiments and 
is 0.946 ± 0.022 (0.991 ± 0.023) times the prediction using 
the Huber+Mueller (ILL+Vogel) fissile antineutrino model. In 
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addition, the measured and predicted spectra are discrepant 
with a significance of ~4er in the 4-6 (5-7) MeV region of 
the IBD prompt (antineutrino) energy spectrum. Investigation 
of the discrepancy strongly disfavors explanations involving 
detector response or an unknown background. A reactor an¬ 
tineutrino spectrum was extracted from the measurement at 
Daya Bay, enabling model-independent predictions of reactor 
antineutrino spectra. 
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